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ABSTRACT 

We have combined the large SNIa database of the Canada- France-Hawaii Telescope Supernova 
Legacy Survey and catalogs of galaxies with photometric redshifts, VLA 1.4 GHz radio sources, and 
Spitzer infrared sources. We present eight SNela in early-type host galaxies which have counterparts 
in the radio and infrared source catalogs. We find the SNIa rate in subsets of radio and infrared 
early-type galaxies is ~l-5 times the rate in all early-type galaxies, and that any enhancement is 
always < 2a. Rates in these subsets are consistent with predictions of the two component "A+B" 
SN la rate model. Since infrared properties of radio SN la hosts indicate dust obscured star formation, 
we incorporate infrared star formation rates into the "A+B" model. We also show the properties of 
SNe la in radio and infrared galaxies suggest the hosts contain dust and support a continuum of delay 
time distributions for SNela, although other delay time distributions cannot be ruled out based on 
our data. 

Subject headings: supernovae: general — galaxies: radio, starburst 



1. INTRODUCTION 

Type la supernovae (SNela) are currently the 
best cosmological tool to investigate the dark en- 
ergy driving the accelerated expansion of the uni- 
verse (Riess et al. 1998; Perlmutter et al. 1999; 
Astier et al. 2006; Wood-Vasey et al. 2007; 

Riess et al. 2007). Though generally accepted as 
thermonuclear explosions of carbon-oxygen white dwarfs 
which have accreted sufficient mass from their compan- 
ion, two competitive models exist for the companion 
scenario: the single-degenerate (main sequence or red 
giant companion) and double-degenerate (another white 
dwarf). Since the time delay between star formation and 
SNela explosion is scenario-dependent, rates in hosts 
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of varying star formation rate (SFR) help constrain 
this. In fact, SNIa rates are known to be correlated 
with host morphology, B-K color, and by inference also 
with star formation (Mannucci et al. 2005). As such, a 
galaxy's SN la rate is commonly expressed as the sum of 
a "delayed" component from old stellar populations and 
a "prompt" component from young stellar populations. 
These components are parametrized as "A" and "B", 
proportional to a galaxy's mass and SFR respectively 
(Scannapieco and Bildsten 2005; Sullivan et al. 2006a). 
This "A+B" two-component model can be matched 
with progenitor populations that have distributions of 
delay times as described by Mannucci et al. (2006), and 
Pritchet et al. (2008). 

The C99 catalog matched with Very Large Array 
(VLA) 1.4GHz radio data revealed the specific SN la rate 
may be enhanced by 2-7 times in radio-loud (Li.4qhz > 
10 29 ergs s _1 Hz -1 ) over radio-quiet (Li.4GHz < 4 x 
10 27 ergs s _1 Hz -1 ) early-type galaxies (Delia Valle et 
al. 2005, hereafter DV05). Although winds from ac- 
tive galactic nuclei (AGN) could increase accretion rates 
of the interstellar medium (ISM) onto white dwarfs, 
a process thought to trigger classical nova eruptions 
(Livio et al. 2002; Madrid et al. 2007), DV05 reject this 
explanation for the enhanced specific SN la rate in radio- 
loud galaxies. They find galaxy interactions are most 
likely to cause the radio emission and supply the neces- 
sary SNIa progenitors via stellar capture during dwarf 
accretion or star formation induced by major mergers 
(Delia Valle et al. 2005). 

The existence of bright infrared counterparts for ra- 
dio galaxies has been well documented and often at- 
tributed to dust obscured star formation coeval with the 
AGN (Magliocchetti et al. 2008; Mainieri et al. 2008). 
By matching a hypothetical galaxy interaction and AGN 
activity timeline with a recurring star formation model 
of ten 10 8 year long episodes each separated by 10 9 
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years, Mannucci et al. (2006) find the enhanced spe- 
cific SN la rate in radio-loud early-type galaxies is best 
fit by a bimodal delay-time distribution (DTD) in which 
the "prompt" (B) delay time is constrained to just < 10 8 
years. They suggest this implies two physical populations 
of SN la progenitors, but note a broad single-population 
DTD could not be ruled out. As this controversial im- 
plication relics on a rate enhancement found with 21 
photometrically identified SNela from DV05, we look to 
confirm this in the large database of spectroscopically 
typed SNe la from the Canada- France Hawaii Telescope 
(CFHT) Supernova Legacy Survey (SNLS). 

We use photometric redshift galaxy catalogs 
(Ilbert et al. 2006) and 1.4 GHz radio catalogs to 
calculate the SNIa rate in radio-loud early- type 
galaxies, and also infrared source catalogs to look for 
obscured star formation in the radio- loud SNIa host 
galaxies. Dust extinction in starbursts hinders SN 
detection: only one SNIa has been detected, and the 
SNIa rate remains unconstrained in starburst galaxies 
(Mannucci et al. 2003; Mannucci et al. 2007). We 
use the infrared catalogs to calculate the SNIa rate 
in bright and luminous infrared galaxies, which are 
known to experience bursts of star formation up to 
100 — 1000 M Q y _1 . We also statistically compare the 
SNIa rates in radio and infrared galaxies to predictions 
of the two-component "A+B" model, and perform an 
identical analysis with the low redshift C99 catalog. 

In § 2 we describe the SNLS and C99 SN la databases 
and galaxy catalogs, the overlapping radio and infrared 
source catalogs, and present the radio and infrared prop- 
erties of SNLS and C99 SNIa early- type host galaxies. 
We derive the SNIa rate in radio and infrared galaxies 
in § 3, compare to expectations of the two-component 
"A+B" model in § 4, and reanalyze with relaxed data 
constraints in § 5. In § 6 we juxtapose the properties 
of the SNela in radio and infrared host galaxies with 
known correlations between SNIa properties and stellar 
populations. Finally, in § 7 we discuss the implications 
of our results regarding current SN la science and future 
surveys, and conclude in § 8. 

2. OBSERVATIONS 

Here we describe the combination and processing of 
private and public catalogs of supernovae, galaxies, and 
radio and infrared sources. 

2.1. CFHT Supernova Legacy Survey Catalog 

From 2003 to 2008, CFHTLS Deep component imaged 
four 1 deg 2 fields (D1-D4) every three to four nights 
(during dark/grey time, when visible) in four MegaCam 
filters (gMji~MiiM, zm) to a depth iu — 25; this was 
accompanied by a strong spectroscopic campaign to fol- 
low up on as many potential SN la candidates as possible 
(Perrett et al. 2009). After final reductions of all data, 
the SNLS 14 will have discovered and identified hundreds 
of SNe la, provide the best direct constraints on the dark 
energy equation of state parameter w (Astier et al. 2006) 
and be a useful compliment to other cosmological tech- 
niques (Komatsu et al. 2009). We use the 91 and 64 
SNela identified prior to 2006 Dec 31 in CFHTLS deep 

14 http://cfht.hawaii.edu/SNLS 



fields 1 and 2 respectively (~ 30% of the final four-field 
total), because suitable VLA radio surveys cover Dl and 
D2 only, and up to this date the SNLS spectroscopic 
analysis and typing is complete. 

2.2. Ilbert et al. (2006) Galaxy Catalog 

The Ilbert et al. (2006) galaxy catalog 15 incorporates 
VIMOS VLT Deep Survey spectroscopic redshifts to cal- 
ibrate their spectral energy distribution (SED) fitting 
routine, resulting in accurate (<ta z /(i+z) — 0.029) photo- 
metric redshifts for galaxies in the four SNLS Deep fields. 
In optimizing the photo-z calculation, the accuracies of 
the SEDs are compromised (Ilbert, private communica- 
tion) and the distribution of SED types is discontinuous. 
To solve this we fit 51 SEDs, interpolated from Coleman 
et al. (1980) and Kinney et al. (1996) templates made by 
Stephen Gwyn at the Canadian Astronomical Data Cen- 
tre 16 , to the catalog galaxies, and estimate galaxy stellar 
masses and star formation rates using fits of this library 
of SEDs to the models of Buzzoni (2005). We correct 
these for systematic offsets (of about a factor of 2) to 
agree with the PEGASE models (Sullivan et al. 2006b). 
Our SED types include E/S0 (ellipticals or early- types), 
Sbc and Scd (spiral galaxies), Irr (irregular galaxies), and 
SB (starbursts, sSFR > 30 x 10" 10 y" 1 ). To ensure cat- 
alog purity we restrict the catalog as recommended by 
Ilbert et al. (2006), and limit to i M < 25. 

2.3. Additional Galaxy and SN Catalog Processing 

We identify SNIa hosts as the nearest Ilbert catalog 
galaxy unless there are two within 2 arcseconds, in which 
case redshift is used to discriminate. In Dl and D2, 22 
and 15 SNe la have no catalog galaxy within 5 arcseconds 
(maximum host offset) and cannot be used; this mainly 
includes SNe la in hosts which do not meet our catalog re- 
strictions (17 and 12), but also a few SNIa whose hosts 
are not in the original galaxy catalgs (5 and 3; many 
are near field edges, or close to foreground stars). For 
SNela with Ilbert catalog hosts, iterative outlier rejec- 
tion is applied to the residual dispersion between host 
photometric and SNIa spectroscopic redshifts for each 
deep field, resulting in photometric redshift uncertain- 
ties of (Ti = 0.032 and 02 = 0.028, rejecting 2 and 6 
SNe la hosts as outliers in Dl and D2. We use SN la with 
spectroscopic redshifts z < 0.6, the limit to which the 
SNLS SN la sample is nearly complete (Neill et al. 2006). 
These restrictions yield 33 and 29 usable SNela in Dl 
and D2 respectively; possible effects of the redshift re- 
striction are considered in § 5. 

For every galaxy we calculate its SNIa observed rate 
per year, Ri a , from the two-component "A+B" model ap- 
plying the time dilation correction to our observed frame 
of reference, as shown in Equation 1: 

A x M + B x SFR 

Ria = \~\~Z ' (1) 

where M is stellar mass (M© ) , SFR is star formation rate 
(M Q / yr), and the A and B values are from Sullivan et 
al. (2006): A = 5.3 ± 1.1 x 10~ 14 SNe y" 1 Mq 1 and 

15 http://terapix.iap.fr 

16 http://www.cadc.hia.nrc.gc.ca/community/CFHTLS- 
SG/docs/cfhtls.html 
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B = 3.9 ± 0.7 x 1CT 4 SNe y -1 (M Q y^ 1 ) -1 . The rate 
Ri a is the number of SNe la expected in a galaxy per 
year of observing. We use it to determine an effective 
control time, C, to account for SNLS survey efficiencies 
such as detection and spectroscopic completeness, and 
the lengths of SNLS observing seasons (approximately 6 
months), for fields Dl and D2 separately (i.e. to have Ci 
and C2) as follows: 



C 



N obs 



la 



(2) 



where N Q b s is the total number of SNLS SNe la observed 
in the field prior to 2006 Dec 31 with z < 0.6, and 53 R-ia 
is the total number of SNe la expected in all z < 0.6 
galaxies per year of observing the field. We find Ci = 
1.580 and C2 = 1.284 where C has units of years. 

2.4. VLA Radio Catalogs 

The VLA-VIRMOS 1.4 GHz Deep survey covers Dl 
to a Si.4GHz flux of 80 /iJy, with a mean rms noise a ~ 
17 ^Jy, and 75% completeness at fluxes Si.4ghz =80- 
180 /uJy (Bondi et al. 2003). Of the 1054 radio sources, 
503 optical counterparts are identified as the closest 
galaxy within 2 arcseconds (maximum VLA positional 
error). Away from the galaxy catalog's masked regions 
around foreground stars our recovery rate is ~ 60%, sim- 
ilar to Ciliegi et al. (2005). The VLA-COSMOS 1.4 
GHz Large Project covers D2 to Si.4GHz ~ 45 ^Jy with 
a mean rms noise of a ~ 15 yuJy (Schinnerer et al. 2006), 
for which we match ~ 50% of radio sources with a galaxy. 

Radio luminosities are derived from photometric red- 
shifts of galaxy counterparts, include a (1 + z) _1 
bandpass correction (Hogg et al. 2002), and are plot- 
ted in Figure 1. A galaxy is radio-loud if Li.4GHz > 
10 29 ergs s -1 Hz -1 , the faint-end limit of the radio 
luminosity function, although exact limits for radio- 
loud, -faint, and -quiet galaxies are not universal 
(Zamfir et al. 2008). Ledlow & Owen (1996) find 14% 
(±2.4%) of elliptical galaxies with absolute R magnitude 
Mr < —20.5 are radio-loud, and amongst C99 elliptical 
galaxies DV05 find 12% (±2%) are radio-loud. In our 
sample of SED type E/S0 galaxies with z < 0.6 we find 
4-8% of galaxies with My < —20 are radio-loud, sug- 
gesting we can identify half of the radio-loud population 
amongst optically bright galaxies due to our higher min- 
imum radio flux limit. 

The apparent bimodality of VLA-COSMOS (D2) radio 
luminosities in Figure 1 is a result of combining the in- 
tegrated fluxes of resolved and unresolved sources. The 
flux of the bimodality valley corresponds to the lower 
limit of integrated fluxes for resolved sources, Si.4GHz ~ 
0.08 mJy, as shown in the right-hand plot of Figure 17 
from Schinnerer et al. (2006). Coincidcntally, it also cor- 
responds to the VLA-VIRMOS (Dl) lower limit. This 
may indicate VLA-VIRMOS and VLA-COSMOS sam- 
ple slightly different radio source populations, with more 
faint, unresolved sources in D2. We consider any influ- 
ence of this on our results in § 5. 



The 



2.5. Spitzer Infrared Catalogs 
Spitzer Wide-area Infrared Extragalactic 



(SWIRE) survey covers Dl to a flux of Ss^m ~ 6.6 [iJy 
(where subscript on S denotes wavelength), and to 



S 24p m ~ 300 ^Jy (Lonsdale et al. 2003). SWIRE ob- 
jects with S3. ejum > 200 yuJy and stellarity > 0.9 are most 
likely stars or QSOs and rejected from the catalog 17 . 
The Spitzer Cosmic Evolution Survey (S-COSMOS) 
covers D2 to S3.6^ m ~ 0.8 ^(Jy and S24 Mm ~ 300 /uJy 
after required aperture corrections (Sanders et al. 2007). 
Objects flagged as likely compromised by nearby bright 
sources are rejected 18 . Both catalogs are available 
via the NASA Infrared Space Archive 19 . We reject 
foreground objects and QSOs, identify optical coun- 
terparts as in § 2.4, and limit to z < 0.6 as in § 2.4. 
The fraction of galaxies infrared counterparts in all 
four IRAC bands and MIPS is ~ 2 — 3%; these are 
plotted on infrared color-color diagrams in Figure 3 
which also show AGN boundaries derived from Spitzer 
First Look Survey data of comparable flux-density 
limits (S3. ejum ~ 7 /xJy and 824pm ~ 300 jitJy) and 
redshifts (z < 0.7) (Lacy et al. 2004). Sajina et al. 
(2005) show the region on the IRAC color-color plot 
(left) is actually appropriate for AGN of redshifts — 2. 
It should be noted that not all galaxies in these regions 
are AGN-dominatcd at IR wavelengths. 

We convert S3 6^m to K-band stellar mass via Balogh et 
al. (2007), and 

S24pm to bolometric infrared luminosity 
Lir and infrared star formation rate SFRir with publicly 
available templates and codes 20 (Chary & Elbaz 2001). 
Plots of Lir versus redshift in Figure 2 show our sample 
of Luminous Infrared Galaxies (LIRG, Lir > 1O 11 L ) 
is complete to z ~ 0.5. The conversion of Lir to 
SFRir is appropriate for starburst-dominated infrared 
emission, but IR color-color diagrams in Figure 3 show 
some IR sources may be AGN-dominated; uncertain- 
ties introduced by applying this conversion to potentially 
AGN-dominated sources are discussed in § 4.2. Figure 
4 shows the optical and infrared masses agree very well, 
but SFR op t is 3-10 times lower than SFRir. Having 
sampled the brightest infrared galaxies which are most 
likely to have obscured star formation, this factor is and 
should be higher than the usual SFRir ~ 2 x SFR opt 
(Takeuchi et al. 2005). 

2.6. Properties of SNLS SNIa Host Galaxies 

We find eight SNLS SNIa host galaxies with radio 
and/or infrared counterparts and z < 0.6, the proper- 
ties of which are given in Table 1 and plotted in Fig- 
ures 1 to 4. Image stamps centred on SNIa coordinates 
are shown in Figure 5. Beyond the redshift limit of 
z = 0.6, there are no radio SNIa hosts, but one with an 
infrared counterpart: the Sbc type host of SNIa 04D2ca 
at z = 0.835 is classified as a LIRG. In § 5 we discuss 
the effects of extending our redshift range to z ~ 1.0. 
Of these eight SNIa hosts, five have SED type E/S0 
and will be included in the analysis: four are radio-loud 
(Li.4GHz > 10 29 ergs s _1 Hz -1 ), three are bright infrared 
galaxies (BIRG, Lir > 1O 1O L ), and two are luminous 
infrared galaxies (LIRG, Lir > 10 11 L Q ). 

Galaxies which are bright at mid-infrared wavelengths 
(MIR, for example 24^m) are described as being star- 

17 http://irsa.ipac.caltcch.edu/data/SPITZER/SWIRE/docs/dclivcry_doc_r2. 
page 43 

18 http://irsa.ipac.caltech.edu/data/COSMOS/gator_docs/scosmos_irac_colDi 

19 http://irsa.ipac.caltech.edu 

20 http://david.elbaz3.free.fr/astro_codcs/chary_elbaz.html 
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burst or AGN dominated, meaning they contain dust 
which absorbs high-energy photons from either young 
stars or AGN, and re-emits them in the MIR. It may 
seem the two radio- loud SNIa host galaxies with SEE) 
type E/SO (normally quiescent) and infrared counter- 
parts must be AGN dominated, but there are several rea- 
sons why this is not the case. First, the radio-loud E/SO 
hosts do not lie along the "AGN plume" of the IRAC 
color-color plot (although 04Dljg and 05Dlhn hosts are 
in its "AGN zone"), and are outside the "AGN zone" 
of the IRAC+MIPS color-color plot as shown in Fig- 
ure 3. Second, AGN-dominated infrared emission orig- 
inates from a hot dusty torus tens of parsecs across 
(Tristram et al. 2009) and would be unresolved, but the 
infrared sources associated with these two galaxies have 
are nagged as extended. Third, a visual inspection of 
radio and infrared SNIa hosts in um and GALEX (D2 
only) images 21 show all but the late-type host of 06D2ff 
are quite faint to invisible in the ultraviolet, which can 
be indicative of dust obscured star formation. Fourth, 
not all starbursts appear morphologically irregular (e.g. 
Figure 4.68 of Binney & Merrifield (1998), in which a 
profoundly disturbed galaxy appears as an elliptical in 
a shallow exposure). We therefore accept the conversion 
of Ltr into SFRir for these radio-loud E/SO SNIa host 
galaxies, as discussed in § 2.5. 

2.7. Cappellaro et al. (1999) Catalog 

The Cappellaro et al. (1999) catalog is a combina- 
tion of visual and photographic searches from Cappel- 
laro et al. (1993) and Evans et al. (1989), provided 
to us by E. Cappellaro. SN data includes name, type 
(classified photometrically), and host name; galaxy cat- 
alog includes Hubble type, recession velocity, B-band 
luminosity and control time. For radio counterparts 
we use the same methods, matching criteria, flux lim- 
its, and radio source catalogs as DV05: the NRAO 
VLA Sky Survey (Condon et al. 1998), the Parkes-MIT- 
NRAO survey (Griffith & Wright 1993), and the GB6 
survey (Gregory et al. 1996). We convert radio flux to 
luminosity as in § 2.4. Table 2 shows we do not recover 
the same radio luminosities but do classify the same hosts 
as radio-loud, except for that of SN 1983J which DV05 
classify as "borderline" radio-loud and thus add 0.5 SN la 
to the number of radio-loud hosts. We find a lower radio 
power and would not do this. 

Galaxy classifications from the NASA Extragalactic 
Database 22 (NED) for the early-type SNIa host sam- 
ple of DV05 (their Table 2) are reproduced in Table 
2, showing some are peculiar (possibly merging) and/or 
have emission lines. We derive mass and star forma- 
tion rate from morphology and optical magnitude as 
in § 2.2, and K-band stellar masses from the 2 Mi- 
cron All-Sky Survey (Skrutskie et al. 2006) as in § 2.5. 
We compile infrared data from two Infrared Science 
Archive 23 catalogs, the IRAS Faint Source Catalog 
(Moshir et al. 1992) and the IRAS Cataloged Galaxies 
and Quasars (Fullmer & Lonsdale 1989), and convert to 
infrared luminosity as described by Sanders & Mirabel 

21 http://www.cadc.hia.nrc.gc.ca/community/CFHTLS- 
SG/docs/cfhtls.html 

22 http: / /ncdwww. ipac.caltech.edu 

23 http://irsa.ipac.caltech.edu 



(1996). We also calculate the expected number of SNIa 
for each galaxy from the two-component "A+B" model 
as described in 2.3, and use the given given galaxy control 
times (the time during which an exploding SN would be 
detectable) which also account for SNe detection biases 
(Cappellaro et al. 1997). Optical properties have an un- 
certainty of ~ 40-60% from conversions between galaxy 
type conventions, magnitude bands, and mass and SFR 
models. 

The properties of early-type C99 SN la hosts presented 
in Table 2 suggest they may not all be simply quiescent 
ellipticals, perhaps containing more SNIa progenitors (of 
the 'prompt' or 'B' component variety) than previously 
expected. For example, NGC 1275 is a LIRG and likely 
contains dust-obscured star formation and a young stel- 
lar component despite being morphologically early-type. 

3. SNIA RATES IN RADIO-LOUD AND 
INFRARED-BRIGHT GALAXIES 

In this section we present SNLS and C99 SNIa rates 
in a variety of early- type host galaxy subsets: radio- 
loud (L 1.4 ghz > 10 29 ergs s -1 Hz -1 ), radio-loudest 
(Li.4GHz > 10 30 ergs s _1 Hz -1 ), bright infrared galaxies 
or BIRG (Lir, > 1O 1O L.0), and luminous infrared galaxies 
or LIRG (Lir > 10 n L o ). The number of SNLS SNIa ob- 
served in Dl-2 prior to 2006 Dec 21 for each subset, N Q b s , 
is shown in Table 3. Also shown (in brackets) is N G b s /C, 
where C is the effective control time from Equation 2 
and accounts for SNLS detection efficiencies and observ- 
ing seasons. Thus, N b s /C is the number of SNIa which 
explode per year of SNLS observations. Also quoted are 
Poisson uncertainties on N b s /C at the 0.84 (la) confi- 
dence level from Gehrels (1986). The amount of mass in 
galaxies in each subset is given in Table 4. We calculate 
the specific SNIa rate in each galaxy subset, sSNRi a , via: 

° SNR " - mwhr ,3) 

The (1 + z) -1 in the denominator is a time dilation 
factor for each galaxy, which converts the units for 
sSNRi a (SNy -1 Mg 1 ) from observed to rest-frame time. 
Specific SN la rates are given in Table 5 in the commonly 
used "supernovae unit" SNuM, which is equal to one su- 
pernova per century per 1O 1O M0. Uncertainties are prop- 
agated from those in Table 3 but do not include the un- 
certainty in galaxy masses. Ratios of SNIa rates in radio 
and IR subsets over all E/SO galaxies are given in Table 
6. These tables also contain equivalent values for the C99 
SN la and galaxy catalogs. We use radio luminosities de- 
rived in § 2.7 to classify field galaxies, but initially use 
radio luminosities from DV05 to classify SNIa hosts in 
order to make a direct comparison with their results. We 
revise our results with our own luminosities and discuss 
in § 3.1 and 4.1. 

3.1. Radio-Loud Galaxies 

With SNLS data, we find the specific SNIa rate in 
radio-loud early-type galaxies to be O.IO2I0 J5 3 SNuM, 
within la of the C99 result in Table 5. Whereas DV05 
quote a ~ 4 times SN la rate enhancement in radio-loud 
over radio-quiet early-type galaxies, VLA radio catalog 
incompleteness means we can identify only half of the 



SN la Rates in Radio Loud Galaxies 



5 



radio-loud population (§ 2.4) and cannot isolate a radio- 
quiet set for the SNLS (see Figure 1). In lieu of compar- 
ing with a radio-quiet set we consider the ratio between 
radio- loud and all early- types as a lower limit, and find it 
consistent with no enhancement at the la and 2a levels 
for the SNLS and C99 catalogs respectively, as shown in 
Table 6. 

Can we correct our radio incompleteness, and does it 
affect our results? While we could use radio-loud galaxy 
mass functions to estimate the mass missing from our 
radio-loud sample, we could not know how many SNIa 
this mass hosted, so cannot correct the SNIa rates for 
radio incompleteness. Instead, we simply use the limited 
mass in radio-loud galaxies we do find, and the limited 
sample of SNIa hosted by it, as a representative sam- 
ple. Any SNIa rate enhancement in radio-loud galaxies 
would also appear in this limited sample, assuming it 
is not proportional to Li.4qhz (i-e. AGN power). The 
SNIa rate is not expected to be proportional to Li.4qhz 
because AGN activity is not suspected of directly influ- 
encing the probability of an SN la explosion, as discussed 
in § 1; any evolution in the physical properties of AGN 
is not expected to influence our results either. Our radio 
incompleteness is thus unlikely to affect our results - they 
are a suitable comparison sample to C99. 

As an extra test, we present the specific SNIa 
rates amongst the "radio-loudest" (Li.4gh z > 
10 30 ergs s -1 Hz -1 ) early- type galaxies in Table 5, 
and their ratio over all early- types in Table 6. Unlike 
the radio-loud, the radio-loudest sample from SNLS 
catalogs is complete to z = 0.6. While the SNLS SNIa 
rate in the radio-loudest early-type galaxies shows no 
increase over the rate in all early-types, the C99 catalog 
shows a ~ 4x enhancement which is twice that for the 
radio-loud sample. However, this is within 3a of a null 
result, and additional factors may artificially inflate this 
rate as discussed below. 

Although we do not recover precisely the same radio 
luminosities for C99 early- type hosts as DV05, we have 
initially used their classifications in order to directly com- 
pare with their results. As shown in Table 2, we would 
place NGC 1316 in the radio-loud subset but not the 
radio- loudest, and would not classify NGC 3106 as "bor- 
derline" radio-loud. Also, SN1968A in NGC 1275 is 
listed as type-I, not type-la. Although rare, core-collapse 
supernovae have been observed in early- type galaxies; 
this is often due to misclassification of late-types as early- 
type, especially for galaxies with HI and radio emission 
(Hakobyan et al. 2008; Bazin et al. 2009). NGC 1275 is 
a known star forming early-type galaxy, has the highest 
Lir of all C99 early- type hosts (Table 2), and is pre- 
cisely the sort to potentially host a CCSN. Consider- 
ing these factors, we have instead 5, 8, and 20 SNIa 
hosts in the radio-loudest, radio-loud, and all early-type 
samples, lowering their rates to 0.123±g;g||, 0.0851^30, 
and 0.042t°;°j l ) 9 SNuM respectively. Ratios for the radio- 
loudest and radio-loud categories drop to 2.9^.1% and 
2.0+^0, bringing the null result within 2a and la re- 
spectively. 

3.2. Bright and Luminous Infrared Galaxies 

With SNLS data, we find the specific SNIa rate 
in early-type BIRG and LIRG to be 0.168,±% 2 ^ 7 and 



0.2491^3? SNuM respectively. As presented in Table 5 
the C99 data results agree within la. Again we consider 
the ratios between B /LIRG and all early- types as a lower 
limit on any possible enhancement, as given in Table 6. 
It appears the SN la rate is a factor of 3-5 times higher in 
B/LIRG, but this is not more significant than la (2a for 
C99 BIRG). While we could have isolated all early- type 
non-LIRG for a ratio, the LIRG fraction is small enough 
that this would not be significantly different from the all 
early-types sample. 

As discussed in § 3.1, SN 1968A in early-type LIRG 
NGC 1275 may not have been a type la. Also, SN 1983J 
and SN 1991Q do not have a definitive classification. The 
scenario that these three are not type la drops the rate 
in C99 BIRG early-type galaxies to 0.041±g;g|| SNuM, 
consistent with the SNIa rate in all early- type galaxies. 
We also note that although one SNIa host is optically 
classified as starburst (see § 2.2), it does not have an IR 
counterpart and cannot be included here. 

4. COMPARISON TO THE TWO-COMPONENT "A+B" 
MODEL 

Specific SNIa rates in subsets of radio and infrared 
early- type host galaxies are ~l-5 times the rate in all 
early- type galaxies, and all enhancement have low sig- 
nificance (2a at most). In this section we test whether 
these rates are also consistent with the two-component 
"A+B" model, and whether the potential enhancements 
might simply be due to star formation providing addi- 
tional SNIa progenitors. To do this, we statistically 
compare the observed number of SN la in radio and in- 
frared early-type galaxies to that predicted by the two- 
component "A+B" model. For the SNLS catalog, the 
total number of SN la predicted by "A+B" over the sur- 
vey's duration is Na+b = C x Ri a for an individual 
galaxy, where Ri a and C are explained in § 2.3. For 
the C99 catalog, Ri a is converted to Na+b using given 
observational control times. Where N Q b s is the number 
observed, the Poisson probability of observing x = N Q b s 
given an expected number /i = Na+b is expressed by 
(Bevington & Robinson 2003): 

P P (x;n)= l ??-^. (4) 

When N b s > Na+b , the summed Poisson probability of 
observing x = N Q b s or more is the integral of this from 
x = N obs to x = oo; when N Q bs < Na+b the probability 
of observing x = N Q b s or less is the integral from x = 
to x = Nobs- It assesses whether the observed number of 
SNIa is consistent with the "A+B" model in any given 
galaxy subset, as used for early- type cluster galaxies in 
Graham et al. (2008). Summed probabilities < 0.05 are 
considered statistically significant results. 

We calculate summed Poisson probabilities for the 
same SNLS and C99 radio and infrared galaxy subsets 
described in § 3. Where optical masses and star for- 
mation rates are used in the "A+B" model, we refer to 
the total expected number of SNIa as NA+B, op t and the 
associated summed probability as P op t- We can also de- 
rive Na+b,ir and Pir by using infrared SFR if a galaxy 
has a MIPS counterpart (using optical SFR if not), an 
account for contributions from dust-obscured star for- 
mation to the "prompt" component. When substitut- 
ing SFRir into the "A+B" model in Equation 1, we 
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must alter the "B" value of Sullivan et al. (2006) be- 
cause it was derived from optical star formation rates, 
and typically SFR m « 2 SFR opt (Takeuchi et al. 2005). 
Although Figure 4 shows SFR m ~ 3-10 x SFR opt , this 
is for the small fraction of brightest infrared galaxies in 
which this factor is, in most cases, due to dust-obscured 
star formation. Had SFRir been used for all regular 
galaxies, the derived B value would be Bir ~ B opt /2. 

4.1. Radio-Loud Galaxies 

Summed Poisson probabilities for SNLS radio-loud and 
radio-loudest early-type galaxy samples are > 0.05 as in 
Table 7, indicating no significant deviation between ob- 
servations and "A+B" model predictions. As discussed 
in § 3.1, it is unlikely the SNLS results are affected by 
our radio incompleteness. 

However, results for the C99 samples suggest signifi- 
cant to very significant deviations in the samples of radio- 
loud and radio-loudest early-type galaxies. These results 
are unchanged with the consideration of infrared SFR. 
Could this be due to a ~ 5 x variation between derived B 
values from different surveys (Greggio et al. 2008)? For 
the C99 radio-loud sample, using 5xB yields NA+B,opt = 
5.9 and P op t = 0.11, but 20xB is necessary to raise 
NA+B,o P t ~ 4 and P opt > 0.05 for the C99 radio-loudest 
early- type subset, so a change in B value cannot be the 
whole story. As discussed in § 3.1, the actual number of 
SN la in the radio- loud and -loudest subsets of C99 may 
be 8 and 5 respectively This would increase P op t and 
Pir to > 0.1 for both subsets, which we do not consider 
significant. 

4.2. Bright and Luminous Infrared Galaxies 

Summed Poisson probabilities in Table 7 for SNLS 
BIRG and LIRG early-type samples hints at a possible 
excess of SNela over optical "A+B" model predictions, 
and the inclusion of SFRir does raise Na+b.ir and Pir, 
but these results are not statistically significant. The 
same can be said for C99 BIRG, though as discussed in 
§ 3.1 the number of C99 SNIa in BIRGs may be just 1, 
in good agreement with NA+B,o P t- In general we find the 
number of SNIa in BIRG and LIRG is consistent with 
predictions of the "A+B" model. 

We note in § 2.5 that infrared emission of poten- 
tially AGN-dominated galaxies has been converted to 
SFRir with a starburst-dominated template, thereby ar- 
tificially inflating N exPj iR. Identifying galaxies which 
lay along the "AGN plume" (left plot, plume between 
-0.1 < log(S 5 . 8 /S 3 . 6 ) < 0.5) and in the "AGN Zone" 
(right plot) of Figure 3 as the most likely to be AGN, we 
find this population contributes 20-40% of the additional 
SNIa predicted when SFRi R is included, as in 20-40% 

of (N A +B,IR)-(N A +B,opt)- 

5. ALTERING THE DATA CONSTRAINTS 

To consider whether altering the constraints on SNLS 
data effects our results we try, in turn: extending our 
rcdshift range to z = 1.0, including galaxies with SED 
type Sbc, and different combination methods for our ra- 
dio catalogs. 

Sixteen SNLS SNela early-type host galaxies have 
0.6 < z < 1.0, but due to the flux limits of our radio 
catalogs, none are associated with a radio or infrared 



counterpart (only the Sbc type host of SNLS SNela 
04D2ca at z = 0.835 is detected as a LIRG). We de- 
termine the effective control time C for Dl and D2 in 
redshift ranges 0.6-0.8 (Ci = 0.897 and C 2 = 0.298) and 
0.8-1.0 (Ci = 0.369 and C 2 = 0.134) to account for vary- 
ing SNLS detection efficiencies. The specific SNIa rate 
in all early-type galaxies to z = 1.0 is ~ 0.054 SNuM, not 
significantly different from that for z < 0.6. Extending 
to z = 1.0 adds mass - but not SNIa - to the radio and 
infrared early-type subsets, reducing rates and ratios to 
insignificant levels and increasing summed Poisson prob- 
abilities to P > 0.10. No new information is gained from 
extending to z = 1.0 as radio and IR catalogs are not 
deep enough, and doing so only dilutes our sample. 

There are 34 galaxies with SED type E/S0 to Sbc 
and z < 0.6, 6 radio-loud and 5 BIRG as in Table 1. 
The specific SNIa rate amongst all E/SO-Sbc galaxies 
is 0.0611^12 SNuM, (higher than for all E/S0); in the 
radio- loud set it is 0.101+^047 SNum (similar) and in the 

BIRG set is O.081±2:oli SNuM (lower). No statistically 
significant SN la rate enhancement in radio or infrared 
over all E/SO-Sbc galaxies is observed, and in these sets 
the summed Poisson probability is P 3> 0.05. Including 
Sbc type galaxies adds to our total SNIa in radio and 
infrared hosts, but their naturally larger B component 
makes it harder to tell whether any contribution to the 
B component is being made from star formation associ- 
ated with radio and infrared emission. 

The deeper flux limit of the VLA-COSMOS in D2 may 
result in slightly different radio source populations be- 
tween Dl and D2. To investigate whether this effects 
our results, we first restrict VLA-COSMOS to the same 
flux limit as VLA-VIRMOS, Si. 4G Hz > 0.08 mJy. Only 
2 SNIa in radio-loud early-type galaxies are observed, 
which decreases the specific SNIa rate to ~ 0.055 SNuM, 
increases its ratio over all early- types to ~ 1, and raises 
the summed Poisson probability to P ~ 0.5. Secondly, 
we limit our analysis to the deeper radio catalog in D2 
and are left with one radio-loud E/S0 SNIa host, which 
yields a rate of ~ 0.085 SNuM, a ratio over all E/S0 
~ 1.7, and P = 0.38. In general we conclude that any 
differences in the radio source populations are not major 
contributors to the results of this experiment. 

6. SN IA PROPERTIES 

It is well established that SN la light curve shape and 
peak luminosity are correlated: more slowly-declining 
light curves reach a brighter peak luminosity. Light 
curves are parametrized by Amis, the decrease in 
magnitude over the first 15 days after maximum light 
(Phillips 1993), or by stretch, s, the amount that a tem- 
plate (average) light curve must be "stretched" to fit the 
observations where s > 1 for brighter, slower declining 
SNela (Perlmutter et al. 1997). SNIa peak luminosity 
depends on the mass of 56 Ni synthesized during the 
explosion (Arnett 1982); whether this mass is influenced 
more by the metallicity or age of the white dwarf pro- 
genitor star is under much debate (Timmes et al. 2003; 
Ropke & Hillebrandt 2004; Gallagher et al. 2005; 

Gallagher et al. 2008; Howell et al. 2009). 

SNIa properties are also correlated with qualities of 
the host galaxies, where late-type star-forming galaxies 
host mostly bright, slowly declining SNela, and early- 



SN la Rates in Radio Loud Galaxies 



7 



type elliptical galaxies with little to no star-formation 
host faint, rapidly declining SNela (Hamuy et al. 1996; 
Howell et al. 2001). Gallagher et al. (2005) compiled a 
database of 57 local SNela and found the m ean A mi 5 in 
late-type (Sa to Peculiar) host galaxies was Am 15 ~ 1.1, 
but in early-type (E/SO) hosts Ami5 ~ 1.45. Sullivan et 
al. (2006) show the median stretch of SNLS SNe la in 
late-type star-forming galaxies is greater than in early- 
types. This is supported by Howell et al. (2007) who use 
light curve fitting routine SiFTO (Conley et al. 2008) 
to constrain the stretch distribution of "B" (prompt) 
component SNe la (associated with young stellar pop- 
ulations) to be centred at s — 1.071 with a = 0.063, 
and of "A" (delayed) component SNe la (associated with 
old stellar populations) to be centred at s — 0.945 with 
a = 0.077. 

If SNela in radio-loud early-type hosts are associ- 
ated with young stars, their stretch values should be 
consistent with the "B" component stretch distribution. 
SiFTO stretch values for SNLS SNela in radio galaxies 
are given in Table 1. (SNIa 06D2je was not included in 
the SNLS third-year results, but preliminary fits show its 
stretch is similar to 05Dlhn). Stretches for SN la 04Dljg, 
05Dlhn, and 06D2ck are consistent with the "B" com- 
ponent distribution of Howell et al. (2007). Individually 
05Dlhn and 06D2ck are more consistent with originat- 
ing from the "A" component at the la level, but collec- 
tively these three have a mean stretch of s = 1.052 with 
a =0.057, which overlaps the la uncertainty of the "A" 
component stretch distribution. 

Amongst C99 SNela in radio-loud early-type hosts, 
five of eight have a Ami 5 given in Table 2 of DV05, 
Ami 5 = {1.28,1.33,1.73,1.88,1.30}. Their mean is 
Anii5 = 1.5 with <TA mis = 0.25 which is consis- 
tent with Ami5 ~ 1.45 for SNela in early-type hosts 
(Gallagher et al. 2005). Instead of a numerical compar- 
ison, DV05 find that the distribution of decline rates 
for SNela in radio- loud early- type hosts is intermediate 
between that found for late-type and early-type hosts. 
They suggest this supports a continuum of SN la and host 
galaxies properties, from young progenitors in very active 
galaxies, intermediate age progenitors in active galaxies, 
and old progenitors in passive galaxies. Stretch values for 
the three SNLS SNe la in radio-loud early-type galaxies 
also support this. 

Two SNLS SNe la found in radio-loud late-type hosts, 
05Dlby and 06D2ff, have stretch values consistent with 
the "B" and "A" component distributions respectively 
which is normal since late-types have an underlying old 
stellar population. Two SNLS SNe la in non-radio BIRG 
hosts, 04D2bt (early-type) and 06D2ca (late-type), have 
stretch values consistent with the "A" and "B" compo- 
nent distributions which also agree with their host type. 

In general we find the properties of radio and infrared 
SNIa hosts, and their environments, to resemble dusty 
starburst galaxies (§ 2.6), and expect their SNIa light 
curves to be reddened and extinguished. SNLS analysis 
finds 05Dlby, 05Dlhn, 06D2ff, and 04D2bt are particu- 
larly red, and the first three have large Hubble residu- 
als, indicating this is an unusually red, faint sample of 
SNela. It is interesting that 04Dljg and 06D2ca are not 
particularly red or faint, because these two SNIa have 
the brightest infrared hosts. Perhaps they originated on 



the 'near' side of the galaxy and escaped dust extinction 
(neither have large host offsets, 2 and 0.5 arcseconds re- 
spectively). Also, the host of 06D2ca is SED type Scd 
and likely has unobscured star forming regions. 

As a final note, SNIa 05Dlby is also spectroscopically 
peculiar in that the 6150A Silicon line is very broad, 
indicating fast ejecta velocities. It is similar to SN 2001ay 
(Nugent et al. 2001), which had an atypically slow light 
curve decline given its unremarkable peak magnitude, 
a very luminous near-IR magnitude at ~ 10 days post 
maximum, and an early-type host (Phillips et al. 2003). 
Such peculiar spectra comprise ~ 1% of all SNela, but in 
this small sample it is impossible to comment on whether 
they are more or less common in active galaxies. 

7. DISCUSSION 

With the C99 and SNLS catalogs we recover the SN la 
rate in radio-loud early- type galaxies of DV05, and for 
the first time we calculate SNIa rates in the radio- 
loudest, BIRG, and LIRG subsets of early- type galaxies. 
We find these rates are ~l-5 times the rate in all early- 
type galaxies, at that any potential enhancement is < 2a. 
For the first time, we have incorporated infrared star 
formation rates into the two-component "A+B" model. 
While it does slightly increase the number of SN la pre- 
dicted, more so in B/LIRG than in radio-loud(est) sam- 
ples, ultimately all SNLS results are consistent with the 
two-component "A+B" model. The same cannot ini- 
tially be said for all C99 results, but the observed num- 
ber of C99 SNela may be an upper limit due to SNIa 
and/or host misclassifications. 

The benefit of considering infrared host galaxy prop- 
erties to account for possible dust obscured star forma- 
tion is best shown in Figure 6, which plots the specific 
SNIa rate as a function of mean specific star forma- 
tion rate (sSFR) for each early- type galaxy subset. Blue 
symbols indicate where sSFR was derived from optical 
data only, and red symbols where SFRir was substi- 
tuted for all galaxies with counterparts in the Spitzer IR 
catalogs. Curved lines show the two-component "A+B" 
model from Sullivan et al. (2006). Although error 
bars are large and all rates are within 2a of the opti- 
cal "A+B" model, when plotted with sSFR m the SNLS 
SNIa rates in B/LIRG and radio-loud(cst) subsets do 
align more closely with "A+B". The same can be said 
for C99 B/LIRG, and the apparent discrepancy for the 
C99 radio-loud(est) subset may be due to misclassifica- 
tions discussed in § 3.1. These results are not significant, 
but arc included to demonstrate the utility of such an 
analysis for future SN la rates. 

DV05 suggested the SN la rate enhancement in radio- 
loud early-type galaxies may be due to galaxy mergers 
triggering radio activity and providing extra SNIa pro- 
genitors by induced star formation. Some aspects of the 
SNLS data are consistent with this: the evidence for 
dust-obscured star formation in SNIa radio hosts, the 
agreement between observations and "A+B" model pre- 
dictions in SNIa radio hosts, and the properties of SNIa 
in radio hosts being similar to those associated with in- 
termediate age stellar populations. However, the dis- 
crepancy between observations and "A+B" model pre- 
dictions for C99 radio galaxies remains. In § 3.1 and 
§ 4.1 we discuss possible SNIa and/or host misclassi- 
fications, but might it be the two-component model is 
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simply inappropriate for radio galaxies? It is a linear ap- 
proximation to a relation that is not necessarily linear, 
and uses current SFR as a measure of the SFR when 
the white dwarf progenitor star was born - the parame- 
ter actually related to the number of SN la progenitors. 
In regular spiral galaxies where the SFR remains con- 
stant for 10 9 to 10 10 Gyr this is essentially true, but not 
necessarily for galaxies experiencing interactions and/or 
mergers, and perhaps episodic star formation (discussed 
below), on timescales of ~ 10 8 years (DV05). This is not 
a failure of the "A+B" model, merely a limitation of its 
application to individual galaxies. In large samples this 
effect should approximately "average out" . 

With such a hypothetical episodic SFR, Mannucci et 
al. (2006) show the radio-loud SNIa rate enhancement 
of DV05 is best fit with a bimodal delay-time distribution 
(DTD) for SN la in which ~ 50% of the prompt compo- 
nent explodes in < 10 8 years, and suggest two physically 
distinct populations of SN la progenitors. However, we 
find the SNLS SNIa rate in radio galaxies agrees with 
the "A+B" model, which is consistent with continuous 
DTD's (Pritchet et al. 2008), and that the properties of 
SNe la in radio hosts do not constrain them to very young 
stellar populations (as does DV05). 

As a final note, the suggestion by DV05 that galaxy 
interactions and/or mergers cause both the radio activ- 
ity and enhanced SNIa rate is consistent with the well 
documented correlation between clustered environments 
and radio galaxies (Magliocchetti et al. 2007). Figure 5 
shows image stamps of radio SNIa host environments; 
galaxies within Az < 0.1 (i.e. ~ 3(Ta z /(i+ z )) are marked. 
This figure demonstrates that most radio loud SNIa 
hosts have nearby neighbours, as expected for elliptical 
galaxies. As reported in Graham et al. (2008), the host 
of 05Dlby is located in the outskirts of a galaxy cluster. 
We also find the galaxy density within r < 50kpc of this 
host shows clustering over the background field galaxy 
distribution, as determined by a "significance parame- 
ter" (Graham et al. 2008); the small-scale environments 
of 04Dljg and 06D2ff hosts also show clustering, but not 
more so than the average SNIa environment. In general 
we find the properties of radio SN la hosts and their en- 
vironments are consistent with dusty starburst galaxies 
in non- isolated regions. A full analysis of SNIa in small 
groups and pairs using spectroscopic (i.e. higher accu- 
racy) galaxy redshifts will be presented in an upcoming 
paper. 

8. CONCLUSION 

Based on the SNLS catalog, SNIa rates in radio and 
infrared early- type galaxies are ~l-5 times the rate in all 
early-type galaxies, but any enhancement is < 2a. Rates 
in these subsets are consistent with predictions of the 
two component "A+B" SNIa rate model. The infrared 
properties of SN la radio-loud early-type host galaxies 
suggest the presence of dust obscured star formation, 
and we have for the first time incorporated SFRir in 
the "A+B" model. In the C99 catalog, radio-loudest 
SNIa hosts are consistent with the "A+B" model only 
if some SNe la and/or host galaxies have been misclassi- 
fied. In general we find the properties of SNIa in radio 
and infrared galaxies support a continuum of SN la and 
host galaxies properties, from young progenitors in very 
active galaxies to old progenitors in passive galaxies, as 



did DV05. Also these SNIa are fainter and redder than 
other SNIa, consistent with a dusty environment, and 
one was spectroscopically peculiar. 

To continue the investigation of SNIa rates in differ- 
ent stellar populations, observations in the CFHT Deep 
fields at near- to far-infrared and sub-millimeter would 
be very useful. Radio coverage of the remaining two 
CFHTLS deep fields could double the number of high- 
rcdshift SN la in radio-loud early-type hosts available for 
analysis, but requires > 50 hours per field to achieve 
the depth of VLA-VIRMOS (Dl) and > 250 hours per 
field for VLA-COSMOS depth (D2). Future work in- 
cluding spectroscopic galaxy redshifts will allow an anal- 
ysis of the small scale environments of SN la hosts, and 
SN la rates and properties in small groups and interact- 
ing pairs. 
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Fig. 1.— Li.4 G Hz versus photometric redshift for galaxies in Dl with VLA-VIRMOS (left) and D2 with VLA-COSMOS (right) coun- 
terparts. Filled circles are galaxies of SED type E/SO (black), Sbc (purple), and later types (blue); squares are SNIa hosts. Solid red line 
marks survey Si.4GHz flux detection limit (dotted red for limit of other field), dashed red line for radio-loud limit. 
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Fig. 3. — IR color-color diagrams for z < 0.6 galaxies with IRAC counterparts only (left) and IRAC+MIPS counterparts (right), for Dl 
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Fig. 4. — IR versus optical mass (left) and SFR (right) for z < 0.6 galaxies IRAC+MIPS counterparts, for Dl and D2 combined. Filled 
circles are galaxies of SED type E/SO (black), Sbc (purple), and later types (blue). Triangles are SNIa host galaxies, with squares for radio 
hosts. Solid lines of slope equal to one are there to guide the eye. 




Fig. 5. — Image stamps in rMegai 30x30 arcseconds, centred on SNela for host galaxies detected in radio and/or infrared source catalogs. 
Right to left, top to bottom: 04Dljg, 05Dlby, 05Dlhn, 06D2ck, 06D2ff, 06D2je, 04D2bt, and 06D2ca. White crosses mark SNela, and 
white plus signs mark nearby galaxies within Az < 0.1 (i.e. ~ 3<J^ Z /(l+z) )• 
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Fig. 6. — Specific SNIa rates in SNLS (filled) and C99 (open) early-type galaxies (circles) from Table 5, versus mean sSFR of each subset. 
Left plot shows radio-loud (RL, square) and radio-loudest (RLcst, star) subsets; right plot shows BIRG (inverted triangle) and LIRG 
(triangle) subsets. Symbol color indicates sSFR derived from optical data only (blue), or by incorporating infrared data when available 
(red). Also plotted is the two-component "A+B" model (solid) with uncertainties (dashed) from Sullivan et al. (2006) (blue lines), and 
the theoretical infrared "A+B" model (with Bjr ~ B op t/2) (red lines) as discussed in § 4. 
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TABLE 1 

Properties of SNLS SN Ia and their Host Galaxies. 



SN Ia 


SNIa 


SN Ia 


SN Ia 


Host 


Host 


Host 


Host 


Host Radio 


Host IR 


Host IR 


SNLS ID 


Zspec 


Color 


Stretch 


SED 


M v 


Mass 


sSFR 


Luminosity 


Luminosity 


sSFR 














[10«» M ] 


[10- 10 y" 1 ] 


[10 29 erg s- 1 Hz- 1 ] 


[10 10 L ] 


[10- 10 y" 1 ] 


04Dljg 


0.584 


-0.092 


1.023 ±0.023 


E/S0 


-22.99 


43.5 


0.15 


14.5 


46.47 


1.56 


05Dlby 


0.299 


0.736 


1.040 ±0.061 


Sbc 


-21.17 


5.43 


0.82 


2.85 


8.972 


2.66 


05Dlhn 


0.149 


0.334 


1.049 ± 0.020 


E/S0 


-22.08 


20.2 


0.08 


2.82 


12.95 


1.09 


06D2ck 


0.552 


-0.020 


1.083 ±0.048 


E/S0 


-22.50 


32.1 


0.03 


5.30 






06D2ff 


0.345 


0.299 


0.980 ±0.092 


Sbc 


-21.94 


9.19 


1.17 


4.16 


7.871 


1.84 


06D2je 


0.418 






E/S0 


-19.79 


2.16 


0.25 


1.73 






04D2bt 


0.220 


0.222 


0.988 ±0.062 


E/S0 


-21.40 


10.1 


0.08 


0.00 


2.193 


0.35 


06D2ca 


0.531 


0.047 


1.086 ±0.049 


Scd 


-21.03 


0.81 


6.83 


0.00 


33.56 


24.6 



TABLE 2 







Prof 


'ERTIES OF C99 ELLIP' 


ncAL Host Galaxi 


ES. 






SN Ia Name 


Galaxy 


Hubble 


Radio 


NED 


Optical 


Infrared 


Infrared 


(and Type) a 


Name a 


Type a 


Luminosity a 


Galaxy Type b 


sSFR 


Luminosity 


sSFR 








[10 29 erg s" 1 Hz- 1 ] 




[IO" 10 y- 1 ] 


[10 10 L ] 


[y- 1 ] 


1961H (Ia) 


NGC 4564 


-4.9 


<0.008 (<0.010) 


E6 


0.00 






1968A (I) 


NGC 1275 


-1.6 


720.0 (1462) 


cD; pec; NLRG 


0.26 


17.5 


0.63 


1970J (Ia) 


NGC 7619 


-5.0 


0.540 (0.590) 


E 


0.00 






1972J (I) 


NGC 7634 


-2.0 


<0.040 (<0.050) 


SB0 


0.27 






19801 (Ia) 


NGC 4374 


-4.7 


20.00 (31.18) 


El; LINER; Sy2 


0.00 


0.2 


0.02 


1980N (Ia) 


NGC 1316 


-1.9 


1000. (1.752) 


SAB; LINER 


0.25 


0.7 


0.03 


1981D (Ia) 


NGC 1316 


-1.9 


1000. (1.752) 


SAB; LINER 


0.25 


0.7 


0.03 


1981G (Ia) 


NGC 4874 


-4.0 


24.00 (27.98) 


cD 


0.03 






1983G (Ia) 


NGC 4753 


-2.2 


0.530 (<0.015) 


10 


0.26 


0.6 


0.08 


1983J (?) 


NGC 3106 


-1.9 


0.920 (0.692) 


SO 


0.26 


1.7 


0.13 


1986G (Ia) 


NGC 5128 


-2.2 


18.00 (52.39) 


SO; pec; Sy2 


0.25 


1.6 


0.20 


1990M (Ia) 


NGC 5493 


-2.1 


<0.040 (<0.051) 


SO; pec 


0.25 






1991Q (?) 


NGC 4926A 


-1.6 


0.310 (0.325) 


SO; pec?; Sbrst 


0.24 


6.2 


2.86 


1991bg (Ia) 


NGC 4374 


-4.7 


20.00 (31.18) 


El; LINER; Sy2 


0.00 


0.2 


0.02 


1991bi (Ia) 


NGC 5127 


-4.9 


8.700 (9.843) 


E; pec 


0.00 






1992A (Ia) 


NGC 1380 


-1.9 


<0.020 (<0.018) 


SAO 


0.25 


0.4 


0.06 


1992bo (Ia) 


E352-G57 


-1.5 


<0.160 (<0.175) 


SB; pec 


0.27 






1993C (Ia) 


NGC 2954 


-4.9 


<0.080 (<0.105) 


E 


0.00 






1993ah (Ia) 


E471-G27 


-2.0 


41.00 (44.77) 


SBb 


0.26 






1994D (Ia) 


NGC 4526 


-2.0 


0.042 (0.047) 


SAB 


0.26 


0.7 


0.15 


1996X (Ia) 


NGC 5061 


-5.0 


<0.030 (<0.033) 


E0 


0.00 







a From DV05, except radio powers in brackets from this paper. 

b Galaxy types (column 5) from the NASA Extragalactic Database (NED) where E: elliptical; cD: supergiant in cluster; I: irregular; 
LINER: low-ionization nuclear emission-line region; NLRG: narrow line radio galaxy; pec: peculiar; S: lenticular, A (unbarred), B (barred); 
Sbrst: star burst; Sy: Scyfcrt. 



TABLE 3 

Number of SN Ia in Early- Type Galaxies 21 



Sample 




SNLS b 


C99 


All 


20 


(14.01-1;?) 


21 +5.7 


Radio-loud 


4 


(2.8tfj) 


9 5+ 4 ' 2 
y '°-3.0 


Radio-loudest 


1 


(0-6±g;J) 


o+3.9 
°-2.8 


BIRG 


3 


(2-Otfj) 


4 +3.2 
4 -1.9 


LIRG 


2 


(i-3±?;S) 


i+2.3 
1 -0.8 



a Unccrtaintics arc derived from Poisson errors at the lcr confidence level, as in Gehrels (1986). 

b Brackcts contain number of SN Ia corrected for SNLS detection efficiencies to be the number which exploded per year of SNLS observations, 
as in § 2.3. 
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TABLE 4 
Mass in Early- Type Galaxies 21 



Sample 


SNLS 


C99 


All 


38440.4 (4337) 


23315.0 (2079) 


Radio-loud 


3879.3 (165) 


7724.5 (314) 


Radio-loudest 


1638.8 (46) 


2825.5 (92) 


BIRG 


1650.3 (179) 


1805.7 (166) 


LIRG 


747.0 (73) 


469.2 (22) 



a Brackcts contain number of galaxies in sample. 



TABLE 5 
SN Ia Rates in Early- Tyi 



•e Galaxies 3 



Sample 



SNLS 
[SNuM b1 



+0.(118 

-0.014 
-0.103 
'-0.057 
-+0.186 
^ -0.051 

0.168±°;?£ 

-0.469 
-0.191 



C99 
[SNuM b1 



+injT2 

0.010 
+0.045 
-0.032 
0.097 
0.068 
+0.131 
-0.079 
0.488 
0.177 



All 

Radio-loud 
Radio-loudest 
BIRG 
LIRG 



0.051 

0.102; 

0.055 



0.249+ 



0.044 
0.101 
0.197 
0.166 
0.213 



a Uncertainties are derived from Poisson errors at the la confidence level, as in Gehrels (1986). 
b [SNuM] = SN (100 yr)- 1 (10 10 M©)" 1 



TABLE 6 

Ratios of SN Ia Rates in Early- Type Galaxies 3 



Sample 


SNLS 


C99 


Radio-loud 


z - u -1.3 


2.3+H 


Radio-loudest 


1 i+5.3 

1 - 1 -i.o 


4 5+ 4 -' 


BIRG 


O Q+6.9 

°-°-2A 


3.7+^1 


LIRG 


49 +14.2 


4 8 + 15 '- 3 



Uncertainties are derived from Poisson errors at the la confidence level, as in Gehrels (1986). 



TABLE 7 

Statistical Comparison to "A+B" 



SN Ia Rate Model 3 



Sample 




Nobs 


NA+B,opt 


Popt 


Na+b,ir 


Pir 


Radio-loud 


SNLS 


4 


2.16 


0.17 


2.38 


0.22 




C99 


9.5 


4.84 


0.04 


5.10 


0.05 


Radio-loudest 


SNLS 


1 


0.90 


0.59 


1.00 


0.73 




C99 


8 


2.35 


0.003 b 


2.39 


0.003 b 


BIRG 


SNLS 


3 


1.03 


0.09 


1.52 


0.20 




C99 


4 


1.39 


0.05 


1.77 


0.10 


LIRG 


SNLS 


2 


0.44 


0.07 


0.83 


0.20 




C99 


1.0 


0.26 


0.23 


0.43 


0.35 



a Columns defined in § 4. 

b Further considerations raise P > 0.1 in : 



4.1. 



